Fast collective motions of backbone in transmembrane α helices are critical to water transfer of aquaporin.

Fast collective motions are widely present in biomolecules, but their functional relevance remains unclear. Herein, we reveal that fast collective motions of backbone are critical to the water transfer of aquaporin Z (AqpZ) by using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and molecular dynamics (MD) simulations. A total of 212 residue site-specific dipolar order parameters and 158 15N spin relaxation rates of the backbone are measured by combining the 13C- and 1H-detected multidimensional ssNMR spectra. Analysis of these experimental data by theoretic models suggests that the small-amplitude (~10°) collective motions of the transmembrane α helices on the nanosecond-to-microsecond timescales are dominant for the dynamics of AqpZ. The MD simulations demonstrate that these collective motions are critical to the water transfer efficiency of AqpZ by facilitating the opening of the channel and accelerating the water-residue hydrogen bonds renewing in the selectivity filter region.

Mechanistic Insight into the Mechanical Unfolding of the Integral Membrane Diacylglycerol Kinase.

The essential forces stabilizing membrane proteins and governing their folding and unfolding are difficult to decipher. Single-molecule atomic force spectroscopy mechanically unfolds individual membrane proteins and quantifies their dynamics and energetics. However, it remains challenging to structurally assign unfolding intermediates precisely and to deduce dominant interactions between specific residues that facilitate either the localized stabilization of these intermediates or the global assembly of membrane proteins. Here, we performed force spectroscopy experiments and multiscale molecular dynamics simulations to study the unfolding pathway of diacylglycerol kinase (DGK), a small trimeric multispan transmembrane enzyme. The remarkable agreement between experiments and simulations allowed precise structural assignment and interaction analysis of unfolding intermediates, bypassing existing limitations on structural mapping, and thus provided mechanistic explanations for the formation of these states. DGK unfolding was found to proceed with structural segments varying in size that do not correlate with its secondary structure. We identified intermolecular side-chain packing interactions as one of the major contributions to the stability of unfolding intermediates. Mutagenesis creating packing defects induced a dramatic decrease in the mechano-stability of corresponding intermediates and also in the thermo-stability of DGK trimer, in good agreement with predictions from simulations. Hence, the molecular determinants of the mechano- and thermo-stability of a membrane protein can be identified at residue resolution. The accurate structural assignment established and microscopic mechanism revealed in this work may substantially expand the scope of single-molecule studies of membrane proteins.

Discovery of Novel Anti-Resistance AR Antagonists Guided by Funnel Metadynamics Simulation.

Androgen receptor (AR) antagonists are widely used for the treatment of prostate cancer (PCa), but their therapeutic efficacy is usually compromised by the rapid emergence of drug resistance. However, the lack of the detailed interaction between AR and its antagonists poses a major obstacle to the design of novel AR antagonists. Here, funnel metadynamics is employed to elucidate the inherent regulation mechanisms of three AR antagonists (hydroxyflutamide, enzalutamide, and darolutamide) on AR. For the first time it is observed that the binding of antagonists significantly disturbed the C-terminus of AR helix-11, thereby disrupting the specific internal hydrophobic contacts of AR-LBD and correspondingly the communication between AR ligand binding pocket (AR-LBP), activation function 2 (AF2), and binding function 3 (BF3). The subsequent bioassays verified the necessity of the hydrophobic contacts for AR function. Furthermore, it is found that darolutamide, a newly approved AR antagonist capable of fighting almost all reported drug resistant AR mutants, can induce antagonistic binding structure. Subsequently, docking-based virtual screening toward the dominant binding conformation of AR for darolutamide is conducted, and three novel AR antagonists with favorable binding affinity and strong capability to combat drug resistance are identified by in vitro bioassays. This work provides a novel rational strategy for the development of anti-resistant AR antagonists.

Jellyfish-inspired smart tetraphenylethene lipids with unique AIE fluorescence, thermal response, and cell membrane interaction.

Smart lipids with fluorescence emission, thermal response, and polyethylene glycolation (PEGylation) functions can be highly valuable for formulation, image-traceable delivery, and targeted release of payloads. Herein, a series of jellyfish-shaped amphiphiles with a tetraphenylethene (TPE) core and four symmetrical amphiphilic side chains were conveniently synthesized and systematically investigated as smart lipids. Compared with regular amphiphilic TPE lipids and phospholipids, the unprecedented jellyfish-shaped molecular geometry was found to enable a series of valuable capabilities, including sensitive and responsive aggregation-induced emission of fluorescence (AIE FL) and real-time FL monitoring of drug uptake. Furthermore, the jellyfish-shaped geometry facilitated the concentration-dependent aggregation from unimolecular micelles at low concentrations to "side-by-side" spherical aggregates at high concentrations and a unique mode of AIE. In addition, the size and the arrangement of the amphiphilic side chains were found to dominate the aggregate stability, cell uptake, and thus the cytotoxicity of the amphiphiles. This study has unprecedentedly developed versatile smart TPE lipids with precise structures, and unique physicochemical and biological properties while the peculiar structure-property relationship may shed new light on the design and application of AIE fluorophores and functional lipids in biomedicine and materials science.

Effect of Host Cholesterol on the Membrane Dynamics of Outer Membrane Lipids of Mycobacteria.

The ability of Mycobacterium tuberculosis to remain dormant after primary infection represents the prime cause of new TB cases throughout the world. Hence, diagnosis and treatment of individuals hosting dormant mycobacterium is one of the crucial strategies to be adopted for the prevention of Tuberculosis. Among many strategies unleashed by the latent bacterium, one of them is scavenging host cholesterol for carbon source. Cholesterol modifies lipid membranes over many scales and here, its effect on mycobacterial membrane biophysics and the subsequent effect on partitioning of antibiotics into cholesterol- enriched mycobacterial membranes was investigated. Our research showed that cholesterol alters the phase state behavior of mycobacterial outer membrane lipids by enhancing the overall membrane order at the headgroup and acyl chain region and is integrated into both ordered and disordered domains/phases, with a preference for the latter. Exogenous cholesterol further alters the drug partitioning behavior of structurally different drugs, pointing to a larger clinical potential of using more hydrophobic medications to target dormant bacteria.

On the Dynamic Mechanism of Long-Flexible Fatty Acid Binding to Fatty Acid Binding Protein: Resolving the Long-Standing Debate.

Fatty acids (FAs) are one of the essential energy sources for physiological processes, and they play a vital role in regulating immune and inflammatory responses, promoting cell differentiation and apoptosis, and inhibiting tumor growth. These functions are carried out by FA binding proteins (FABPs) that recognize and transport FAs. Although the crystal structure of the FA-FABPs complex has long been characterized, the mechanism behind FA binding and dissociation from FABP remains unclear. This study employed conventional MD simulations and enhanced sampling technologies to investigate the atomic-scale complexes of heart fatty acid binding proteins and stearic acid (SA). The results revealed two primary pathways for the binding or dissociation of the flexible long-chain ligand, with the orientation of the SA carboxyl head during dissociation determining the chosen path. Conformational changes in the portal region of FABP during the ligand binding/unbinding were found to be trivial, and the overturn of the ″cap″ or the unfolding of the α2 helix was not required. This study resolves the long-standing debate on the binding mechanism of SA with the long-flexible tail to FABP, which significantly improves the understanding of the transport mechanism of FABPs and the development of related therapeutic agents.

Fusion Landscape of Mycobacterial Envelope-Derived Lipid Vesicles with Intact Bacteria Dictates High Intracellular Drug Retention.

Membrane vesicles are critical regulators of pathogenic diseases. In tubercular infections, the use of mycobacteria derived vesicles as delivery vehicles to overcome drug resistance and complex treatment regimens has never been attempted. Here, we first address how these vesicles interact with their target cells, especially via membrane fusion. Membrane fusion between alike mycobacterial outer and inner membrane layer-derived lipid vesicles is shown to be driven by the structural, geometrical, and biophysical attributes of constituent lipids. The increased fusion of outer-membrane-derived vesicles with intact bacteria ensures enhanced intracellular drug levels and is presented as a "natural" antitubercular drug delivery vehicle.

The Conformational Transitions and Dynamics of Burkholderia cepacia Lipase Regulated by Water-Oil Interfaces.

Structural dynamics and conformational transitions are crucial for the activities of enzymes. As one of the most widely used industrial biocatalysts, lipase could be activated by the water-oil interfaces. The interface activations were believed to be dominated by the close-to-open transitions of the lid subdomains. However, the detailed mechanism and the roles of structure transitions are still under debate. In this study, the dynamic structures and conformational transitions of Burkholderia cepacia lipase (LipA) were investigated by combining all-atom molecular dynamics simulations, enhanced sampling simulation, and spectrophotometric assay experiments. The conformational transitions between the lid-open and lid-closed states of LipA in aqueous solution are directly observed by the computational simulation methods. The interactions between the hydrophobic residues on the two lid-subdomains are the driven forces for the LipA closing. Meanwhile, the hydrophobic environment provided by the oil interfaces would separate the interactions between the lid-subdomains and promote the structure opening of LipA. Moreover, our studies demonstrate the opening of the lids structure is insufficient to initiate the interfacial activation, providing explanations for the inability of interfacial activation of many lipases with lid structures.

Molecular Insights on the Conformational Transitions and Activity Regulation of the c-Met Kinase Induced by Ligand Binding.

The tyrosine-protein kinase Met (c-Met) is an important signaling molecule involved in cellular growth and division. The dysregulation of c-Met may induce many fatal diseases, including non-small cell lung cancer, gastrointestinal cancers, hepatocellular carcinoma, etc. The activation of the c-Met kinase is dominant by the structure and dynamics of many important functional motifs, which are regulated by adenosine triphosphate (ATP) binding. c-Met inhibitors bind to the ATP-binding site or the allosteric pocket to compete with ATP molecules or alter the conformation of the function-related domains. Nevertheless, the mechanisms of ligand binding to c-Met are still unclear, especially the regulation of the functional motifs by different inhibitors. These greatly impede the development of novel drugs to overcome the drug tolerance to the currently marketed c-Met inhibitors. In this study, we used enhanced sampling technology to study the binding and regulation of two specific c-Met inhibitors. The results show that the two ligands adopt different binding processes even though with similar binding affinity. More importantly, our results uncovered different protein conformational features and the correlated motions of functional motifs regulated by the inhibitors, providing the structural basis for the functional suppression of the protein kinases.

Mutually Exclusive Interactions of Rifabutin with Spatially Distinct Mycobacterial Cell Envelope Membrane Layers Offer Insights into Membrane-Centric Therapy of Infectious Diseases.

The mycobacterial cell envelope has spatially resolved inner and outer membrane layers with distinct compositions and membrane properties. However, the functional implication and relevance of this organization remain unknown. Using membrane biophysics and molecular simulations, we reveal a varied interaction profile of these layers with antibiotic Rifabutin, underlined by the structural and chemical makeup of the constituent lipids. The mycobacterial inner membrane displayed the highest partitioning of Rifabutin, which was located exclusively in the lipid head group/interfacial region. In contrast, the drug exhibited specific interaction sites in the head group/interfacial and hydrophobic acyl regions within the outer membrane. Altogether, we show that the design of membrane-active agents that selectively disrupt the mycobacterial outer membrane structure can increase drug uptake and enhance intracellular drug concentrations. Exploiting the mycobacterium-specific membrane-drug interaction profiles, chemotypes consisting of outer membrane-disruptive agents and antitubercular drugs can offer new opportunities for combinational tuberculosis (TB) therapy.

The Conformational Transition Pathways and Hidden Intermediates in DFG-Flip Process of c-Met Kinase Revealed by Metadynamics Simulations.

Protein kinases intrinsically translate their conformations between active and inactive states, which is key to their enzymatic activities. The conformational flipping of the three-residue conservative motif, Asp-Phe-Gly (DFG), is crucial for many kinases' biological functions. Obtaining a detailed demonstration of the DFG flipping process and its corresponding dynamical and thermodynamical features could broaden our understanding of kinases' conformation-activity relationship. In this study, we employed metadynamics simulation, a widely used enhanced sampling technique, to analyze the conformational transition pathways of the DFG flipping for the c-Met kinase. The corresponding free energy landscape suggested two distinct transition pathways between the "DFG-in" and "DFG-out" states of the DFG-flip from c-Met. On the basis of the orientation direction of the F1223 residue, we correspondingly named the two pathways the "DFG-up" path, featuring forming a commonly discovered "DFG-up" transition state, and the "DFG-down" path, a unique transition pathway with F1223 rotating along the opposite direction away from the hydrophobic cavity. The free energies along the two pathways were then calculated using the Path Collective Variable (PCV) metadynamics simulation. The simulation results showed that, though having similar free energy barriers, the free energy cuve for the DFG-down path suggested a two-step conformational transition mechanism, while that for the DFG-up path showed the one-step transition feature. The c-Met DFG flipping mechanism and the new intermediate state discovered in this work could provide a deeper understanding of the conformation-activity relationship for c-Met and, possibly, reveal a new conformational state as the drug target for c-Met and other similar kinases.

Self-Assembled Oligopeptide (FK)4 as a Chiral Alignment Medium for the Anisotropic NMR Analysis of Organic Compounds.

Anisotropic NMR parameters have been proven to be powerful for the structural elucidation of organic molecules. Herein, we present an alignment medium based on the self-assembled (FK)4 oligopeptide, showing excellent properties in measurements of anisotropic NMR parameters in both D2O and CD3OD. The preparation of the (FK)4-based alignment medium is simple and rapid. The low viscosity of the anisotropic phase makes it easy to be transferred to the NMR tube. The alignment of the oligopeptide is fast, stable, and homogeneous, with weak background signals, permitting the acquirement of high-quality NMR spectra. The performance of this alignment medium in residual dipolar coupling measurements and diastereomer discriminations is demonstrated by analyzing several different analytes. The enantiodiscrimination property of the (FK)4 oligopeptide is revealed by the difference of residual chemical shift anisotropy of the two enantiomers in the 1D 13C spectrum, granting its potential use for the quantification and identification of enantiomers of small molecules.

Lipid Clustering in Mycobacterial Cell Envelope Layers Governs Spatially Resolved Solvation Dynamics.

The mycobacterial cell envelope acts as a multilayered barrier to drugs. However, the role of lipid composition in the properties of different mycobacterial membranes, otherwise dictating their interactions with drugs, is poorly understood. In this study, we found that hydration states, solvation relaxation kinetics, rotational lipid mobility, and lateral lipid diffusion differed between inner and outer mycobacterial membranes. Molecular modeling showed that lipid clustering patterns governed membrane dynamics in the different layers of the cell envelope. By regulating membrane properties, lipid composition and structure modulated water abundance and interactions with lipid head groups. These findings can help deepen our understanding of the physical chemistry underlying membrane structure and function, as well as the interaction of mycobacterial membranes with drugs and host membranes.

The auto-inhibition mechanism of transcription factor Ets-1 induced by phosphorylation on the intrinsically disordered region.

As the most abundant post-translation modifications (PTMs), the phosphorylation usually occurred on the intrinsically disordered regions (IDRs). The regulation on the structures and interactions of IDRs induced by phosphorylation is critical to the function performing. The eukaryotic transcription factor 1 (Ets-1) is a member of transcription factor family, which participates in many important biological processes. The DNA-binding ability of Ets-1 is auto-inhibited by a disordered serine-rich region (SRR) on the Ets-1. The inhibition ability of SRR is greatly enhanced by the phosphorylation of the serine on the SRR. Nevertheless, the molecular mechanisms of the phosphorylation regulation on the structure and activity of Ets-1 are still unclear and under debates. By using both of the molecular simulations and biochemical experiments, we studied the molecule mechanism of phosphorylation regulation on the auto-inhibition of the Ets-1. The reasons of stabilization of Ets-1 core by phosphorylation on SRR region were elucidated. More important, the free energy landscapes (FEL) show that both of the steric hindrance and allosteric regulation are responsible for the DNA-binding inhibitory induced by phosphorylation, but the steric effects contribute greater than the allosteric regulation. The phosphorylation not only enhances the electrostatic interactions to facilitate the steric impedance, but also promotes the formation of hydrophobic residue clusters, which provide major driven force for the allosteric regulation. The structural basis of auto-inhibition of Ets-1 induced by the phosphorylation revealed in this study would great help the developing of inhibitor for the cancer therapy.

Metadynamics Simulations to Study the Structural Ensembles and Binding Processes of Intrinsically Disordered Proteins.

The structures of intrinsically disordered proteins (IDPs) are highly dynamic. It is hard to characterize the structures of these proteins experimentally. Molecular dynamics (MD) simulation is a powerful tool in the understanding of protein dynamic structures and function. This chapter describes the application of metadynamics-based enhanced sampling methods in the study of phosphorylation regulation on the structure of kinase-inducible domains (KID). The structural properties of free pKID and KID were obtained by parallel tempering metadynamics combined with well-tempered ensemble (PTMetaD WTE) method, and the binding free energy surfaces of pKID/KID and KIX were characterized by bias-exchanged metadynamics (BE-MetaD) simulations.

Phosphorylation Regulation on the Homo-Dimeric Binding of Transactive Response DNA-Binding Protein.

The dimerization of transactive response DNA-binding protein of 43 kDa (TDP-43) is crucial for the RNA metabolism, and the higher-order aggregation of TDP-43 would induce several neurodegenerative diseases. The dimerization and aggregation of TDP-43 are regulated by the phosphorylation on its N-terminal domain (NTD). Understanding the regulation mechanism of TDP-43 NTD dimerization is crucial for the preventing of harmful aggregation and the associated diseases. In this study, the dimerization processes of wild-type (WT), phosphorylated S48 (pS48), and phosphomimic S48E mutation (S48E) of TDP-43 NTD are characterized by the enhanced sampling technology. Our results show that the phosphorylation not only shift the conformation population of bound and unbound state of TDP-43 NTD, but also would regulate the dimerization processes, including increase the binding free-energy barrier. The phosphomimic mutation would also shift the conformational space of TDP-43 NTD dimer to the unbound structures; however, the thermodynamic and kinetic properties of the dimerization processes between the phosphorylated and phosphomimic mutant systems are distinct, which reminds us to carefully study the phosphorylation regulation by using the phosphomimic mutations.

Molecular View on the Dissociation Pathways and Transactivation Regulation Mechanism of Nonsteroidal GR Ligands.

As a major drug target for anti-inflammatory therapy, the glucocorticoid receptor (GR) regulates a wide range of physiological processes through transactivation (TA) or transrepression. GR TA is involved in many adverse effects of GR-targeting drugs, and therefore, the discovery of novel GR ligands with lower TA activity and longer residence time is quite urgent. Undoubtedly, understanding the ligand dissociation mechanisms and the structural basis of the TA regulation is crucial for the development of novel GR-targeting drugs. Here, we used random accelerated molecular dynamics (RAMD) and funnel metadynamics (FM) simulations to explore the dissociation mechanisms of 5 classic glucocorticoids and 6 nonsteroidal GR ligands. Multiple ligand dissociation pathways were discovered. The classic glucocorticoids exhibit a strong preference for Path I, and most nonsteroidal ligands tend to dissociate along mixed pathways. We also find that the distinct unbinding preferences for AZD2906 and AZD9567, two representative nonsteroidal ligands with similar scaffolds but different TA activities, are primarily determined by their different polar interactions with the surrounding residues. Notably, the binding of AZD9567 poses a substantial impact on the conformation of the GR homodimer interface, which provides a valuable clue to understand the mechanisms of the TA-related side effects induced by the adjustments of the homodimerization process. These findings are critical for the structure-based rational design of novel GR ligands with more potent anti-inflammatory potency and reduced side effects.

Discovery of Novel GR Ligands toward Druggable GR Antagonist Conformations Identified by MD Simulations and Markov State Model Analysis.

Binding of different ligands to glucocorticoid receptor (GR) may induce different conformational changes and even trigger completely opposite biological functions. To understand the allosteric communication within the GR ligand binding domain, the folding pathway of helix 12 (H12) induced by the binding of the agonist dexamethasone (DEX), antagonist RU486, and modulator AZD9567 are explored by molecular dynamics simulations and Markov state model analysis. The ligands can regulate the volume of the activation function-2 through the residues Phe737 and Gln738. Without ligand or with agonist binding, H12 swings from inward to outward to visit different folding positions. However, the binding of RU486 or AZD9567 perturbs the structural state, and the passive antagonist state appears more stable. Structure-based virtual screening and in vitro bioassays are used to discover novel GR ligands that bias the conformation equilibria toward the passive antagonist state. HP-19 exhibits the best anti-inflammatory activity (IC50 = 0.041 ± 0.011 µm) in nuclear factor-kappa B signaling pathway, which is comparable to that of DEX. HP-19 also does not induce adverse effect-related transactivation functions of GR. The novel ligands discovered here may serve as promising starting points for the development of GR modulators.

Discovery of a Novel Androgen Receptor Antagonist Manifesting Evidence to Disrupt the Dimerization of the Ligand-Binding Domain via Attenuating the Hydrogen-Bonding Network Between the Two Monomers.

Androgen receptor (AR) has proved to be a vital drug target for treating prostate cancer. Here, we reported the discovery of a novel AR antagonist 92 targeting the AR ligand-binding pocket, but distinct from the marketed drug enzalutamide (Enz), 92 demonstrated inhibition on the AR ligand-binding domain (LBD) dimerization, which is a novel mechanism reported for the first time. First, a novel hit (26, IC50 = 5.57 µM) was identified through virtual screening based on a theoretical AR LBD dimer bound with the Enz model. Then, guided by molecular modeling, 92 was discovered with 32.7-fold improved AR antagonistic activity (IC50 = 0.17 µM). Besides showing high bioactivity and safety, 92 can inhibit AR nuclear translocation. Furthermore, 92 inhibited the formation of the AR LBD dimer, possibly through attenuating the hydrogen-bonding network between the two monomers. This interesting finding would pave the way for the discovery of a new class of AR antagonists.

Regulation Mechanism for the Binding between the SARS-CoV-2 Spike Protein and Host Angiotensin-Converting Enzyme II.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is mainly mediated through the interaction between the spike protein (S-pro) of the virus and the host angiotensin-converting enzyme II (ACE2). The attachment of heparan sulfate (HS) to S-pro is necessary for its binding to ACE2. In this study, the binding process of the receptor-binding domain (RBD) of S-pro to ACE2 was explored by enhanced sampling simulations. The free-energy landscape was characterized to elucidate the binding mechanism of S-pro to ACE2 with and without HS fragment DP4. We found that the stability of the T470-F490 loop and the hydrophobic interactions contributed from F486/Y489 in the T470-F490 loop of S-pro are quite crucial for the binding, which is enhanced by the presence of DP4. Our study provides valuable insights for rational drug design to prevent the invasion of SARS-CoV-2.